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Abstract. We discuss some of the main open issues related to the light-up and evolution of 
the first accreting sources powering high redshift luminous quasars. We discuss the perspec- 
tives of future deep X-ray surveys with the International X-ray Observatory and possible 
synergies with the Wide Field X-ray Telescope. 
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1. Introduction 

The "dark ages" of the Universe ended when 
the UV radiation from the first objects reion- 
ized the intergalactic medium. This epoch be- 
gan approximately at 6 < z < 14 and so far has 
remained essentially unexplored. Investigating 
the end of the dark ages is extremely important 
to understand structure formation, since dur- 
ing this phase the first proto-galaxies and the 
first seed black holes, that would later grow 
into luminous quasars, formed. Observations 
of the most distant (z> 6) quasars (QSOs) thus 
probe the very early growth of supermassive 
black holes (SMBHs) in the centers of massive 
galaxies. 

The high luminosities and broad line 
widths of the most distant QSOs r equire 
BH masses above 10 9 M o dFan etal.1 120011: 
IWillott et al.l 120031) . The formation of such 
systems within the first 700-800 Myr of cos- 



mic time is a challenge for theoretical models 
ri.e. IVolonteri & Rees 1 120051: iLi et al. I l2007t 
Begel manT 20101) . A number of possibilities 
have been proposed for the origin of the 
seed BHs (from Pop III stars to more mas- 
sive objects resulting from the direct col- 
lapse of molecular clouds), the accretion rate 
(Eddington limited or super Eddington) and 
the merging rate of dark halos in the early 
Universe. Even though suitable combinations 
of the above parameters may explain the pres- 
ence of massive BHs at z ~ 6, the processes 
responsible for their assembly and light-up are 
largely unknown. The search for and the study 
of the first QSOs is a key scientific goal of fu- 
ture X-ray missions and in particular of the 
Wide Field X-ray Telescope (WFXT; Rosati et 
al. 2010, this book) and the International X-ray 
Observatory (IXO; White et al. 2010). 
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Fig. 1. The AGN number density, as determined by recent X-ray surveys, in two 2-10 keV luminosity 
ranges. Left panel: logLx=43 -44.5; right panel: logLx=44.5-47. Green circles and squares at z < 3 are from 
the XMM-COSMOS survey dBrusaetal.ll2010h . Black circles and squares and open red circles at 2 < z < 
3 are from various surveys as described in Aird et al. (2010). Filled red circles at z > 3 in the right panel are 
from Brusa et al. (2009). The best fit results of the various determinations of the X-ray luminosity function 
from Aird et al (2010; black solid line), Ebrero et al. (2009; pink dashed curves), Silverman et al. (2008; 
red dashed curves) are reported in both panels. The green dot-dashed curves represent the expectations of 
the XRB synthesis model of Gilli et al. (2007). 



2. Supermassive Black Holes in the 
early Universe 

Most of the SMBH accretion luminosity is 
emitted in the optical/UV and X-ray bands. As 
a consequence, optical and X-ray surveys have 
played a major role in the discovery of high- 
redshift quasars. 

2.1. Optical surveys 

A major step forward in the study of the high 
redshift Universe was achieved thanks to the 
Sloan Digital Sky Survey (SDSS). In a se- 
ries of papers, Fan and collaborators searched 
for z > 5.7 QSOs in different releases of the 
SDSS over an area of about 6600 deg 2 with 
a color selection (/-band drop out) technique 
(see Fan et al. 2001 for a detailed descrip- 
tion). Their final sample includes 19 objects 
with 5.7 < z < 6.4 and Mi 450 < -26. Using 
the same technique, the search for z - 6 QSOs 
was recently extended towards lower luminosi- 
ties (M1450 ~ -24) in the deep SDSS stripe 
by Jiang et al. (2008,2009), allowing the dis- 



covery of 12 additional objects. Thanks to the 
Canada France High Redshift Quasar Survey 
(CFHQS), which covers about 500 deg 2 down 
to fainter magnitudes, and applying a similar 
selection criterion, 1 9 additional QSOs at z > 
5.7 ha ve been found dWillott et al. ll2007Ll2009l 
|2010|) . The number of optically selected QSOs 
revealed up to z ~ 6 is large enough to de- 
termine their luminosity function. By combin- 
ing the CFHQS data with t he luminous SDSS 
QSOs, IWillott et alTI (1201 Ol) assembled a sam- 
ple of 40 QSOs to estimate the 5.74 < z < 6.42 
luminosity function over the range -24.7 < 
M]4so < -27.5. Although the covariance be- 
tween the slope and M* 45Q luminosity prevents 
strong constraints to be placed on either pa- 
rameter, there is evidence for the detection of 
a break at M* 45() =* -25.1. As far as luminous 
QSOs are concerned, an exponential decline (oc 
j q-o.43z) m jjjgjj. S p ace density fromz ~ 3 to z ~ 
6 has been clearly observed. Although SDSS 
has pushed the redshift record well above z ~ 
6 with a statistics good enough to estimate the 
cosmological space density of high-z QSOs, 
it is crucial to point out that SDSS QSOs are 
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among the brightest and most extreme sources 
in the early Universe (in terms of luminosities 
and BH masses), and thus probably not repre- 
sentative of the QSO population at those red- 
shifts. 



2.2. X-ray surveys 

The Chandra and XMM sensitivities to faint 
and hard X-ray sources has opened up a new 
era in the observations of the high redshift 
Universe. The systematic study of high-z (> 



4) optically selected QSOs was pioneered b 



z 



Brandt and coworkers (i.e. | Ka sni et al. ]|200 
IVignali et al. 120031; Ishemmer et al. 120061) : see 
Brandt et al. 2004 for a review) with pointed 
snapshot observations, and further expanded 
by cross correlating the SPSS catalogues 
with the XMM archive dYoung et al. I 2009) 
and the Champ project dGreen et al. I [2009). 
Multifrequency observations point towards the 
absence of evolution in the X-ray spectral 
index, absorption, metallicity, emission-line 
strengths and dust properties over a broad 
redshift range suggesting that the physical 
mechanisms powering luminous QSOs are in- 
sensitive to the significant c hanges on larger 
scales that occur at z - 0-6 dMaiorinoll2006t 
IShemmer et al. Il2006h . Some departures from 
a self-similar evolution have emerged thanks 
to Spitzer near infrare d observations of a small 
sample of z - 6 QSOs djiang et al. Il2010h . Two 
out of the 21 z —6 QSOs observed by Spitzer 
do not show any detectable emission from hot 
dust. There is also evidence that hot dust free 
QSOs have the smallest BH masses in the sam- 
ple (2-3 x 10 8 Mq) and are accreting close to 
their Eddington limit, suggesting an early evo- 
lutionary stage. 

A major improvement in the study of the 
cosmological evolution of the X-ray emission 
of moderate to high-z QSOs has been obtained 
in the last few years thanks to the large num- 
ber of deep and medium-deep Chandra and 
XMM surveys. Large samples of X-ray se- 
lected objects were built, and by now the X- 
ray luminosity function is sampled up to z ~ 
3^1, down to X-ray luminosities of the or- 
der of a few xlO 44 erg s _1 and to z ~2-3 
at lower luminosities (see Fig. 1; Brus a et aTl 



| 2009t lEbrero et al. 112009b lYencho et al. Il2009t 
Uird et al. 11201 OtlFontanot et al. 112007^ 

At higher redshifts, present X-ray obser- 
vations are highly incomplete, being strongly 
limited by the survey areas. The redshift 
records for X-ray selected QSOs are reached 
in the COSMOS survey (z = 5.41; Civano 
et al 2010 ApJ in prep.) and in the CLASXS 
survey (z = 5.40; Steffen et al. 2004), with 
only one spectroscopically identified object per 
survey. Deep and ultra-deep Chandra surveys 
in the CDFS and CDFN would be sensitive 
to lower X-ray luminosities, but the lack of 
area coverage is much more severe. Also in 
this case only one z > 5 spectroscopically con- 
firmed AGN is foun d in the CDFN at z = 5.19, 
dBareer et al. I l2005l) . and none in the CDFS 
(ILuoetal.ll2010h . As a result, the space den- 
sity of X-ray selected QSOs is still uncon- 
strained atz > 5. 



3. Beyond z ~ 6 

The volume and surface density of the first 
QSOs can be estimated following two differ- 
ent approaches. A theoretical one, by building 
up semi analytical models (SAMs) which in- 
clude a description as accurate as possible of 
the complex physical processes thought to be 
at work and key model parameters such as the 
seed BH mass, the accretion rate and the peak 
in the density field fluctuations above which 
they can start to form. Alternatively, an obser- 
vational route may be followed. The best fit lu- 
minosity function and evolution, as measured 
by present X-ray surveys or by optical surveys 
can be extrapolated to high redshifts with some 
educated guesses on the optical to X-ray lumi- 
nosity ratio. 

Fig. 1 shows the most updated compila- 
tion of the X-ray luminosity function (XLF) 
parameterized with a Luminosity Dependent 
Density Evolution phenomenological model 
(LDDE; Silverman et al. 2008; Ebrero et al. 
2009) or a Luminosity and Density Evolution 
model (LADE; Aird et al. 2010). The predic- 
tions based on the AGN syn thesis model for 
the XRB of iGilli et al.l d2007l) are also shown, 
obtained by including an exponential cut-off 
in the XLF evolution at z=2.7, which is a good 
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Fig. 2. Left panel: predicted AGN number counts at z > 6 for various semi-analytic models of BH/galaxy 
formation (solid lines): Salvaterra et al. (2007, green), Rhook & Haehnelt (2008, blue) and Marulli et al. 
(2009, red). Dotted curves show extrapolations of the luminosity functions adopted in the Gilli et al. (2007) 
model for the synthesis of the XRB (see text for details). The vertical orange line represents the IXO 
confusion limit (reached in about 2 Ms). Right panel: the sky coverage of a possible "survey" program (see 
text for details). 



description of the observed space density of lu- 
minous QSOs in th e XMM-COSMOS survey 
dBrusa et al. Il2009h . 

A detailed description of the best fit pa- 
rameters of each XLF is beyond the pur- 
poses of this paper. Here we would like to 
stress that the extrapolations beyond the red- 
shift range where the XLF was computed 
may differ by more than one order of magni- 
tude. Moreover, the faint end of the XLF is 
poorly determined, even at relatively low red- 
shifts (left panel of Fig. 1), making the pre- 
dictions even more uncertain. The space den- 
sity of z > 6 QSOs predicte d by a few recent 
SAMs dMarulli et al. 1 120091: ISalvaterra et al. I 
l2007tlRhook & Haehnelt 1120081) are compared 
with the extrapolation of the observed XLF in 
the left panel of Fig. 2. 

The lower dotted curve corresponds to the 
"decline" model (cfr. Fig. 1) which is found to 
provide a good fit to the observed space density 
of z ~ 3-5 QSOs in the Chandra-COSMOS 
survey (Civano et al. 2010, in prep.). The upper 
dotted curve corresponds to an ad-hoc param- 
eterization of the XLF high-redshift evolution 
named "maxXLF". The space density of low 
luminosity (logLx < 44) AGN is kept constant 
at z > 4. This model maximizes the predicted 



number of higher redshift QSOs, being at the 
same time in agreement with the observational 
results at lower z, hence the name "maxXLF" 
(see Fig. 1 in Gilli et al. 2010, this volume). 
Different choices for the evolution of the lumi- 
nosity function (cfr. Fig. 1) would correspond 
to space densities in between the "decline" and 
the "maxXLF" predictions which should then 
be considered as a conservative and an opti- 
mistic estimate, respectively. 



4. Breaking through the first typical 
AGN 

The bulk of the population of high-z QSOs 
representative of the first accreting objects in 
the Universe is likely to be characterized, on 
average, by relatively low X-ray luminosi- 
ties and BH masses. Large area optical and 
near-infrared surveys such as those foreseen 
with future ground-based facilities and space 
observatories like panstarrs, vista, lsst and 
jdem/euclid will discover a large number of 
high-z QSOs. The above surveys will be bi- 
ased against dust obscured QSOs and faint ac- 
tive nuclei for which the host galaxy starlight 
cannot be neglected. 
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Fig. 3. Residuals wrt a single power law fit to IXO simulations. Left panel: a 100 ks simulation, including 
background, of a z = 6 QSO. The line EW of ~ 40 eV observed frame (280 eV rest frame) is clearly seen at 
~ 0.9 keV. Right panel: a 1 Ms simulation, including background, of a Compton Thick AGN at z = 7. The 
characteristic strong iron line (EW ~ 1.2 keV rest-frame) is easily recognizable. 



X-ray observations can probe much lower 
luminosities and high absorption column den- 
sities. Moreover, the X-ray emission is the tell- 
tale of an accreting object. A blind search for 
X-ray selected z ~ 5-6 QSOs would require to 
survey several tens of square degrees to a depth 
comparable to that reached in the Chandra 
deep surveys (about 0.2 square degrees for a 
total observing time of 6 Ms) and is far beyond 
the capabilities of current X-ray facilities. 

Thanks to its large effective area at 1 keV 
(~ 3 m 2 ) and good spatial resolution (5 arcsec 
HPD), IXO will perform deep surveys reach- 
ing limiting fluxes below 10~ 17 erg cirT 2 s _1 
in the 0.5-2 keV energy range (Fig. 2, left 
panel). The final sensitivity will depend on 
the actual number of faint sources which sets 
the confusion limit. In order to harvest a sig- 
nificant number of moderate luminosity, ob- 
scured AGN at z > 6 and to cope with the 
present large uncertainties on their predicted 
space density, a typical multi-cone observing 
strategy is mandatory. 

A possible observational strategy for a 
multi-cone IXO survey with the Wide Field 
Imager (WFI) is shown in the right panel of 
Fig. 2. The total observing time invested in 
cosmological surveys would be of the order 
of 18 Ms broadly split in deep (3 pointings, 2 
Ms each) medium (20 pointings, 300 ks each) 
and shallow (100 pointings, 60 ks each) cor- 
responding to about 10 months of observing 



time assuming an efficiency of 70%, though 
it should be noted that shallower observations 
will be performed anyway for a wide range of 
scientific investigations. 

The expected number of z > 4 QSOs will 
range from a few hundreds to more than a thou- 
sands, while at z > 6 the predicted number of 
sources could vary beween a dozen up to a few 
hundreds, reflecting the uncertainties reported 
in Fig. 1 and 2. Taking a median value, the to- 
tal number of high-z QSOs will be sufficient 
to constrain the evolution of the faint end of 
the XLF down to L x ~ 5x 10 42 erg s -1 and up 
to z -6. IXO will also be able to detect the very 
first objects at z ~ 10, if they exist, and if their 
X-ray luminosity is larger than approximately 
10 43 erg s _1 . IXO deep surveys need to be per- 
formed over well studied sky areas at longer 
wavelengths to ease the process of source iden- 
tification and redshift determination. The ex- 
cellent IXO angular resolution and sensitivity 
coupled with the capabilities of jwst, alma, e- 
elt and tmt will allow to study the assembly of 
early Black Holes and their host galaxies. 

The large IXO collecting area is extremely 
well suited for detailed spectroscopic studies. 
According to the currently accepted models for 
the joint evolution of SMBH and their host 
galax ies ([Hopkins et al. 2006; La mastra et al. I 
2010), the fraction of obscured and heavily ab- 
sorbed or Compton Thick (Nh > 10 24 cirT 2 ) 
AGN is predicted to increase with redshift 
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(i.e. lMenci et al, l20Q8h . From an observational 
point of view, the issue of an increased frac- 
tion of obscured AGN at high redshifts is 
still debated (see the discussion in Gilli et al. 
2010). Irrespective of what the true fraction is, 
the study of heavily obscured AGN at high-z 
would bear important information on the for- 
mation and early growth of the first active 
galaxies. Primordial, gas and dust rich galax- 
ies, should undergo a phase of high obscura- 
tion, gas accretion on the central BH and vig- 
orous starformation, making deep hard X-ray 
observations an ideal tool for their discovery. 
The IXO spectroscopic capabilities are shown 
in the two panels of Fig. 3. The spectrum of a 
luminous (Lx — 3x 10 44 erg s _1 ), mildly ob- 
scured (N H ~ 10 23 cirT 2 ) QSO at z=6 is sim- 
ulated with a 100 ks WFI exposure. The 6.4 
keV iron line EW is assumed to be 280 eV 
rest-frame. The residuals, with respect to an 
absorbed power law fit show a line-like fea- 
ture with an observed frame EW of ~ 40 eV. In 
the right panel, the X-ray emission of a z =7, 
L x ~ 10 43 erg s _1 , CT AGN is modeled with 
a pure reflection spectrum plus a strong (EW ~ 
1 .2 keV rest-frame) iron Ka line. The residuals 
of the fit with the same model, without includ- 
ing line emission, show that the Ka iron line 
can be used to directly measure source redshift 
with an accuracy of Az -0.2. 

5. IXO and WFXT Synergies 

The study of the first SMBH would greatly 
benefit from joint WFXT and IXO surveys. 
The synergies between the two proposed mis- 
sions are obviously clear. The WFXT multi- 
tiered survey strategy is designed to maximize 
the yield in the discovery of high-z AGN. 
Detailed and quantitative estimates are exten- 
sively discussed in Gilli et al. (2010, this vol- 
ume). IXO with its superior sensitivity to faint 
sources would be unique for spectroscopic 
follow-up of the large WFXT samples and 
to discover the faintest AGN at the highest 
redshifts. At the flux of 10~ 15 erg cirT 2 s _1 , 
about 300 AGN at z > 6 are expected in the 
WFXT medium survey assuming a conserva- 
tive "decline" model. Relatively inexpensive 
IXO pointings (100 ks each), would allow us 



to obtain spectra of the quality show in the left 
panel of Fig. 3, making possible to measure 
metal abundances and many other QSO physi- 
cal parameters at very high-z. 
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